Effect of amplification on conductance distribution of a disordered waveguide 
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Introduction of optical gain to a disordered system results in enhanced fluctuations [i*V 2 ) = 
var(g) / (g) 2 ] of dimensionless conductance g, similar to the effect of Anderson localization in passive 
medium. Using numerical simulations we demonstrate that despite of such qualitative similarity, 
the whole distribution of conductance of amplifying random media is drastically different from that 
of passive system with the same value of F(2) ■ 
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Wave transport in quasi- ID (Q1D) is a paradigm of 
mesoscopic physics Universality of the statistics of 
wave transport in Q1D geometry makes it a convenient 
system to study interference effects. Anderson localiza- 
tion phenomenon|2j is an ultimate manifestation of the 
wave interference UQ. Light transport in a Q1D system 
is described by three directly measurable quantities^: 
T a b - the transmission coefficient from an incoming chan- 
nel a to an outgoing channel b; T a = T a b ~ the total 

b 

transmission coefficient from channel a to all outgoing 
channels; g — ^T a (, - conductance. This is unlike elec- 

a,b 

tronic systems where only conductance can be accessed 
experimentally. Landauer formula^ puts the transport 
of electronic Q and electromagnetic waves on the same 
footing. Average value of conductance ((j) e j > 1 
diffusion, g <C 1 - localization), its variance, and more 
generally, its entire distribution has been found to be in- 
dicative of the nature of wave transport through random 
mediumpl • Coherent amplification of light adds a new 
dimension to the study of mesoscopic transport that can 
lead to new phenomena such as random lasing 9] . 

The effect of amplification leads to enhancement of 
nonlocal correlationsllOl and the fluctuations of various 
transport coefficients |l ll 1 1 2j that bears similarity to that 
of localization. Selective amplification of the quasi- modes 
with long lifetime increases their contribution to trans- 
port, leading to "localization by gain" |lfj . lllLll3lll4j . In 
our previous work we found that P(T a b) in active ran- 
dom medium coincides jllj. within numerical accuracy, 
with that of a passive system^] with reduced value of 
dimensionless conductance which parameterizes the dis- 
tribution. The question arises about its relevance to the 
actual value of conductance and its distribution in ran- 
dom medium with gain. 

Scaling theory of localization^] predicts that distri- 
bution of dimensionless conductance, g, depends on a sin- 
gle parameter. The average value (g) is usually chosen 
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for the parameterization!?! @- Q1D geometry (optical 
waveguide) is a particularly convenient system to study 
localization-delocalization transition, which can be real- 
ized with an variation of the length of waveguide. 

Introduction of optical gain to the system increases the 
average conductance of the waveguide. At the same time, 
it is also expected to enhance the conductance fluctua- 
tions. It has been shown that fluctuations of both T ab 
and T a [3 

grow faster than the corresponding averages. 
However, the effect of gain on statistical properties of g 
is not known. 

In Ref. ^3 m or der to obtain correlation functions 
that contain meaningful information about light trans- 
port in amplifying random media, we introduced the con- 
ditional average over all non-lasing realizations (...) — > 
(...) c - Due to fluctuations of lasing threshold such re- 
alizations always exist however improbable. This fact 
makes analytical treatment of the problem difficult. 

In this work, following the same approach [Tol H?| . we 
obtained numerically the full distribution of conductance 
in passive and active disordered waveguides with time 
reversal symmetry We demonstrate that the second 
moment F( 2 ) = var(g)/g 2 = {(g 2 ) — g 2 ) / g 2 increases 
with gain and that the conductance distribution P(g) in 
amplifying system is drastically different from the distri- 
bution of passive system with the same second moment 

In our numerical simulation, we consider 2D waveg- 
uide filled with random medium (Fig. The walls of 
the waveguide are metallic, and circular scattering par- 
ticles are dielectric with refractive index n = 2. We 
use finite difference time domain method to calculate 
the response of our system to pulsed excitation, followed 
by Fourier transformation which gives us the desired 
continuous- wave response^. We have successfully used 
this method to study mesoscopic fluctuations and nonlo- 
cal correlations [lfj, The system we are considering 
is Q1D, and the transition from diffusion to localization 
can be realized by increasing the length L of the random 
medium. The effect of absorption or gain (inside the 
scatterers) is treated by classical Lorenzian model with 
positive or negative conductivity. The advantage of our 
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FIG. 1: Disordered waveguide geometry. 



numerical model is the ability to introduce spatially uni- 
form gain as well as to separate coherent amplification of 
an input signal from spontaneous emission of the active 
medium. In the presence of gain, long after the short ex- 
citation pulse, the electromagnetic field decays with time 
in the non-lasing realizations, while it keeps increasing in 
the lasing ones. We excluded the lasing realizations from 
our statistical ensemble. 

In order to calculate dimensionless conductance, g, a 
short pulse was launched via a series of equidistant point 
sources which uniformly cover the input end of the waveg- 
uide. One calculation was performed for each source in 
the sequence, and the transmission coefficient was ob- 
tained in the series of detector points which uniformly 
cover the output end of the waveguide. The summation 
of T a {, over source and detector points gave us a quan- 
tity which is proportional to dimensionless conductance, 

T a b — OLg. The coefficient of proportionality a is ge- 

ometrical, and is constant for all realizations. It can be 
found as follows 18] . According to scaling theory of local- 
ization [lif , the distribution of conductance in a passive 
system is uniquely determined by its average value. We 
use normalization-free second moment -F( 2 ) of the system 
to find actual value of average conductance g. This can 
be done using results of Rcf. where g and var(g) 
were calculated analytically. Note that once the coeffi- 
cient a is found from passive system (where F( 2 ) and g 
are uniquely related) it can be used in the amplifying 
or absorbing systems with the same arrangement of the 
source and detector points. 

Far into the diffusive regime P{g) has a Gaussian 
shape with the width determined by the universal con- 
ductance fluctuations. In terms of quasi- modes of diffu- 
sive system [l^. addition of weak uniform gain to the sys- 
tem should result in partial compensation of the modes 
decay rates. Significant changes should occur when a 
substantial fraction of the quasi-modes are brought suffi- 
ciently close to their lasing threshold. This is the regime 
which we would like to investigate in our Q1D system, 
that lead to the following choice of system parameters. 
We start with the passive system close to the onset of lo- 
calization with g — 1.53, as determined by the procedure 
outlined above. The waveguide supports N = 20 modes. 

Solid line in Fig. Ufa) shows frequency dependence [2(J 
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FIG. 2: (a) Solid, dash-dotted and dashed lines represent 
g(f) in passive, absorbing and amplifying systems respec- 
tively. Three dashed curves (from bottom to top) correspond 
to the gain coefficient a- g /u gtCr — 0.11, 0.22 and 0.33 (& g ,cr is 
the critical value for average lasing threshold |Tcjp. The ab- 
sorption coefficient for the dash-dotted line o a = — 0.45<7 SiCr .. 
(b) g(f) of a passive system with nearly the same Fn) as the 
upper dashed line of (a). 



of the g in one of random realizations of the passive 
system. As one adds gain (dashed lines) or absorption 
(dash-dot line), the conductance increases or decreases 
respectively. Absorption is known to reduce fluctuations 
of the conductanc el2ll |22| and makes the distribution 
approach Gaussian [23|. The latter fact led to a proposi- 
tion of statistical formulation of localization criterion]^ 
in absorbing systems: parameter g' = (2 /3)var[T a / (T a )] 
takes on the role of the average dimensionless conduc- 
tance g as the localization parameter. It was conjectured 
that localization in absorbing system should be marked 
by the reduction of g' to below unity, similar to g < 1 in 
passive systems. 

Similarly, we can introduce an alternative localization 
parameter based on the fluctuation of conductance. Fig. 
EJa) shows that introduction of gain leads to enhanced 
fluctuations of conductance despite of the increase of its 
average value. We define a new localization parameter 
g" = (2/15)-F(2)j where the numerical factor 2/15 came 
from the leading term in the perturbative expansion of 
variance ofo in the diffusive regime 0,0- Based on data 
from Ref. y| we estimate that in passive systems g" = 1.0 
when g = 0.91. 

In order to make comparative analysis of the full distri- 
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FIG. 3: The inset plots the normalized second moment 
F(2) — var (§)/g 2 as a function of the average conductance g 
in system with different amount of gain (triangles correspond 
to (T g /(Tg,cr = 0.11, 0.22 and 0.45) or absorption (squares 
correspond to (T a bs/&g,cr = —0.11, —0.22 and —0.45). Pas- 
sive system (circle) is shown for reference. In the main plot, 
solid, dash-dotted and dashed lines represent P(g) in passive 
(g — 1.53) or P c (§s) in absorbing ((J a /r g , cr = —0.45) and am- 
plifying (two dashed curves from right to left correspond to 



0.22 and 0.45) systems. Geometrically all systems 



are the same and differ only by the presence/absence of gain 
or absorption. The difference between left-most dashed curve 
and dotted curve is discussed in the text. 



bution of conductance in active medium (c./. Fig. Ufa)) 
with that in passive medium (c.f. Fig. |2Ib)) we first find 
a passive system with the same value of second moment 
Fm) . Then we scale down the conductance of system 
with gain g s — Cg so that (g s ) — g paS sive- Note that 
Fn) = var(g)/g 2 will not be affected by the scaling. 

Figure 13 shows the change of P(g s ) with introduction 
of gain (absorption). From the inset of Fig. 3 we see 
that in amplifying random medium as the optical gain 
coefficient increases, var(g) increases faster than the av- 
erage value of conductance (i.e. normalized second F^ 2 ) 
moment increases with g). This is also evident from the 
shift of the scaled distribution towards smaller conduc- 
tances. The increase of F( 2 ) results in smaller (g s ) by 
virtue of the scaling. Secondly, the shape of the distribu- 
tion also changes, which we will discuss below. The effect 
of absorption is significantly more moderate: the dash- 
dotted (absorption) and the left dashed curve correspond 
to systems with the same absorption/gain coefficients. 

Now we would like to return to the question of defini- 
tion of conditional statistics in amplifying random media 
that we first discussed in Ref. [lfj. Our numerical algo- 
rithm allows calculation of transport coefficients within 
a certain frequency range. This range is chosen narrow 
enough so that there is no appreciable change in physical 



parameters such as transport mean free path There- 
fore, we use the conductance at different frequencies to 
increase the number of entries in our statistical ensembles 
to about ~ 10 5 . When obtaining conditional statistical 
ensemble we dropped the disorder realizations that ex- 
ceeded the lasing threshold. Thus, the conductances at 
the non-lasing frequencies for the lasing configurations 
are dropped, although lasing does not necessarily affect 
the conductances in the entire frequency range obtained 
numerically. Direct calculation of conductance in this 
paper allows us to test the extent of this effect. In Figs. 
13141 we plotted P c {g s ) calculated using two conditional 
ensembles, (i) Disorder configurations (and thus contri- 
bution from all frequencies in that realization) which lead 
to diverging intensity are dropped - dashed lines in the 
figures, (ii) The entire P(g) is plotted first. The obtained 
distribution contained two separate parts, the second one 
at large g is artificial and contains the lasing contribu- 
tions. Its position is determined by the running time of 
our simulation and therefore can be well separated from 
the contribution of nonlasing configurations/frequencies. 
By discarding the contributions above the gap separat- 
ing two maxima in P{g) we obtain the sought conditional 
distribution - dotted curves in Figs. 13141 

Comparison of two conditional distributions shows 
that (ii) leads to slightly greater fluctuations. In case 
of our samples F/ 2 ) j = 0.12 whereas Fi 
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0.13 com- 



pared to F( 2 ) = 0.055 in the passive system. Removal of 
the contribution from lasing configurations at non-lasing 
frequencies reduces fluctuations of conductance and may 
not be well justified. However, it leads to relatively small 
correction as it can be seen from Figs. 13141 

The question arises: how does P c {g s ) compares to 
P(g) of passive system (which exhibits the same fluc- 
tuation, F(2))l Fig. 0|compares the conditional distribu- 
tions of dimensionless conductance in amplifying random 
medium with that in passive medium with the same F( 2 ) 
(g P assive — 0.92) 0. One can see that the distribution 
in the system with gain has markedly narrower peak, but 
significantly more extended tail toward large g s . This is 
in contrast to the intensity distribution|lll| for which we 
found a good fit with the formula^ of passive system 
parameterized with the reduced value of conductance. 
The result presented in Figs. 13141 suggests that strong 
modification of conductance distribution occurs in the 
systems with gain. 

Optical gain is known to enhance the contribution 
from photons which propagate along long paths and 
therefore spend long time inside random medium before 
escaping|26|. These rare long paths should make con- 
tributions to the tails of the conductance distribution of 
passive random medium. This is because longpaths are 
likely to close on themselves and form loops that is 
known to lead to universal conductance fluctuations Q. 
In the language of quasi- modes (eigenmodes of Maxwell's 
equations of a finite-size random medium without gain 
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FIG. 4: Conditional distributions of the scaled dimensionless 
conductance g s (dashed and dotted lines) in random ampli- 
fying medium (a" 9 /<7 9iCr = 0.45) is compared with that in 
passive medium (bold solid line) with nearly the same -F(2)- 
Thin solid and dash-dotted lines (same as in Fig. [3J corre- 
spond to the original passive system and absorbing system 
with absorption length equal to amplification length of the 
system with gain. 

or absorption) , with the introduction of optical gain, the 
quasi-modes with longer lifetime experience more ampli- 
fication, and they generate the large peaks in the conduc- 
tance spectrum [Fig. 2(a)]. Thus they form the large-g s 
tail of P c (g s ). Meanwhile, small-g tail should be greatly 
suppressed. The above conclusions are indeed supported 
by our numerical data, c.f. Fig. Strong reshaping of 
the conductance distribution with an increase of ampli- 
fication also explains why it differs drastically from that 
of passive random medium with the same variance. 

In conclusion, we computed numerically the distribu- 
tion of dimensionless conductance in the systems with 
and without optical gain. The fluctuation of conductance 
is amplified in active random media. We introduced the 
concept of conditional distribution of conductance P c {g s ) 
that omits the contribution of lasing realizations of ran- 
dom structures and scales the conductance for compari- 
son with that of passive system. Unlike conditional inten- 
sity distribution of P^Tai,)^]], the obtained conductance 
distribution can no longer be fitted by that of passive 
medium with reduced value of the average conductance. 
P c {g s ) in an amplifying random medium differs signif- 
icantly from that of the passive system with the same 
value of normalized variance, P( 2 ) = vcir(jg)/0) . It is 
skewed towards large values of g that we attribute to a 
dramatic reshaping of the distribution of quasi-modes. 
In passive localized systems the transport is dominated 
by the modes that show the strongest coupling to the 
outside baths, whereas in amplifying systems the most 



decoupled, confined modes would be preferentially am- 
plified and relocated into large-g tail of the conductance 
distribution. 
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